Feshbach resonances of harmonically trapped atoms 
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Employing a short-range two-channel description we derive an analytic model of atoms in isotropic 
and anisotropic harmonic traps at a Feshbach resonance. On this basis we obtain a new parameter- 
ization of the energy-dependent scattering length which differs from the one previously employed. 
We validate the model by comparison to full numerical calculations for ^Li-*^Rb and explain quan- 
titatively the experimental observation of a resonance shift and trap-induced molecules in exited 
bands. Finally, we analyze the bound state admixture and Landau-Zener transition probabilities. 



In the last decade reams of fascinating experiments 
with ultracold atoms have been carried out with applica- 
tions ranging from studying condensed matter Hamilto- 
nians and new phases of matter to performing quantum 
information processing 'l|. Two key techniques made 
these achievements possible: (i) Atom-atom interaction 
characterized by the s-wave scattering length a can be 
tuned using a magnetic Feshbach resonance (MFR). (ii) 
Atoms can be confined in various geometries such as 
dipole traps, optical lattices, or atomic waveguides [2]. 
The known theory of MFR's successfully describes the 
free scattering process for varying magnetic field B and 
energy _E — > 0. However, for the full understanding and 
precise controllability of confined atoms at an MFR a 
trap- version of this theory is needed which incorporates 
the energy dependence of the scattering process. This 
is especially the case for strong confinement which has 
been lately used to explore confinement-induced reso- 
nances and scattering in mixed dimensions [3|. These 
systems show exciting behavior such as the formation of 
confinement-induced molecules [4]. 

In the following we present an approach to analyt- 
ically describe MFR's of harmonically trapped atoms. 
We obtain the eigenenergy equation, the admixture of 
the resonant molecular bound state A{E,B), and a gen- 
eral parameterization of the energy-dependent scattering 
length a{E,B). We show that the energy dependence 
of a differs significantly from the one previously used 
to describe trapped gases [a, |6| while it confirms the 
functional form of other two-channel models for trapped 
atoms [7[. We derive energy-dependent formulations of 
the resonance width and resonance shift and extend the 
model to anisotropic harmonic traps. The validity of our 
approach is verified by comparing to full multi-channel 
calculations for ®Li-^''Rb in harmonic confinement. 

We demonstrate the usefulness of the new model by ex- 
plaining the experimental observation of a shift of the res- 
onance position of *^Rb in an optical lattice j8| and by cir- 
cumstantiating the observation of confinement-induced 
molecules in exited states by Syassen et al. [9|. We ana- 
lyze the bound state admixture and show that it might 
be responsible for enhanced losses of trapped ^Li far away 
from the resonance [lO|. Finally, Landau-Zener transition 
probabilities are derived from the full energy spectrum. 



We consider the relative-motion s-wave Hamiltonian 

»-2 .2 

of two atoms in harmonic confinement H = — ^ -^ -t- 
^muj'^E? -f VzH + Vi-ai.{R) with reduced mass m, trap 
frequency w, Zeeman and hyperfine energy Vzh, and 
the electron-spin dependent interaction potential Vint (^)- 
Within the two-channel (TC) description one projects 
onto the subspace of open and closed channels with the 
operators P and Q, respectively. We consider the case of 
an elastic collision with one open channel. This results 
in the coupled equations 



(Hp-S)|*p)+W|*q) 


= 0, 


(1) 


(Hq-^)|*q)+W^|*p) 


-0, 


(2) 



with Hp = PHP, Hq = QHQ, W = PHQ, |*p) = P|*), 
I*q) = Ql^)i and E the energy above the threshold of 
the open-channel interaction potential |ll| . Furthermore, 
one assumes that close to the MFR |^q) is simply a mul- 
tiple ^ of a bound eigenstate |${,) with eigenenergy E^. 
We call this closed-channel state "resonant bound state" 
(RES). To first order, the energy Ei, may be expanded 
linearly in the magnetic field i?, i.e. E\j{B) = ^{B~ Bq). 
Be !$£;) the normalized solution of the open channel 
with |*p) = C\^e) then C'^+A'^ = l holds which allows 
us to define a phase tam^nes = A/C attributed to the 
RBS admixture. Introducing |*p) = C|$b) and |*q) = 
y4|${,) into Eq. ^ and multiplying by ($b| gives 



tan Orbs — 



(^g|iy|$fc) 

{E - Eb) 
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For infinite detuning E — Eb the open channel is as- 
sumed to be in some background eigenstate |$bg) of Hp 
with eigenenergy Ehg- Multiplying Eq. ([T]) by |$bg) and 
using Eq. (|3]) yields the eigenenergy equation 



{E - Eb){E - E^,^) = 






(4) 



In order to find simplified expressions for ($bg|I^|'&&), 
(<I'£;|VI^|$b), and {^hgl^E) we assume that the interac- 
tion acts only in some small range R < Rmt "C Oho 
such that for R > i?int the solution \^e} is given by 
$£;(i?) — A,yD^(p), where D^{p) is the parabolic cylin- 
der function, p ~ ^/2R/aho, Oh, 



== y/h/Jmuj), u = 



E/{huj) — 1/2, and A^, is a normalization constant. For 
R < flho one has D^{p) = D^(0) (l - p/[^/2f{E)]) with 
/(^) - r (i - 5I,) /2r (f - J,) [a. Considering the 
logarithmic derivative one obtains the scattering length 
a{E) = ahof{E), which is equivalent to the result in [13|. 
In the spirit of a Taylor expansion we parameterize 
{^E\W\^b) by a linear combination a$_E(0) + ^$^(0). 
That is, one can define a 7_e and an a* such that 



with resonance width AB 



($b|W^|$6> =7£^b(0) 1 



(5) 



Here, 7^; describes the coupling strength to the RBS and 
a* defines the scattering length of the state \^e) when 
it is orthogonal to W^|$f,). Since the orthogonality fixes 
the phase of |$b) within the coupling range the energy- 
dependence of a* should be usually negligible. We find 
that also the variation of je is negligible which can be 
explained by the stability of the nodal structure of |$b) 
for most of the coupling range. Analogous to Eq. ([5|), 

we set ($bg|W|$f,) = 7bg$Bbg(0) (l - ^j where we al- 
low for a different coupling strength 7bg of the uncoupled 
background state to the RBS. We assume a* to have the 
same value as in Eq. ([5]) since it is determined by the re- 
quirement of orthogonality to the constant term W^|$;,)- 

Finally, we set ($bg|*E) = (^bgl^s) neglecting the 
behavior of the wave- functions at i? < i?int. This ap- 
proximation cannot reproduce the exact energies E where 
(i'bgl'&B) = that depend on the nodal structure at 
R < Rint- However, it is applicable in a sufficient range 
around E — _Ebg such that states of any energy can be 
described by choosing an appropriate background state. 

Applying these assumptions to Eq. Q and using the 
properties of {D^\D^') jl4ij one finds the analytic eigenen- 
ergy equation 



27bg7. (/(^) " £) {f(^- 



' a^ohLo f{E) - /(ii;bg) • ^ ' 

In order to determine the interaction dependent scat- 
tering length a{E, B) we demand that it is equal to the 
scattering length ahof{E) of ^^(i?), i.e. that the eigenen- 
ergies are given by the roots of a{E, B) — ahof{E). Since 
^int ^ Oho the trap has no influence on the interaction 
and thus on a{E,B). The value of /(Ebg) in Eq. © 
is in analogy determined by the root of the eigenequa- 
tion for the uncoupled problem abg(£^bg) = aho/(^bg) 
which is closest to E. Here, Obg is the background scat- 
tering length that varies with the energy approximately 



like abg(fc) ^ 
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k RoB with k = 2mE/h , uq the 



zero-energy background scattering length, and ReS > 
the effective range that can be well estimated from the 
van-der-Waals coefficient Cq [15|. Since the interaction 
is trap independent one can use the a; — > limit to set 
f{E\jg) = abg(-E)/aho. Then, rearranging Eq. ([6]) yields 



27bg7i;"iQbf 



(1- -^ 



and 



detuning SB = abgA_B/(abg — a*). The scattering length 
a* ~ a{E — 0, Bq) has an important impact on the be- 
havior of Ai? and 5B. For small Uhg/a* one has AB oc 
a^ and SB « const, while for systems with large Obg the 
value of a* is negligible such that SB = AB ex Obg. 

Let us compare the result to the previously used energy 
dependence of the scattering length given as [a, |6| 



a{E,B) = abg 1 



AB (1 + (/cflbg 



B^Bo + SB+ {kabg^AB ~ E/^iJ 

(8) 
Here, the term (/cabg)^ induces an additional energy- 
dependence. We examined a two-channel model system 
and found no energy dependence connected to (fcabg)^ 
while the behavior described by Eq. ^ could be vali- 
dated [l6j. The absence of the dependence on (/ccbg)^ is 
also supported by other two-channel models in the pres- 
ence of a trapping potential [7] . 

From the formal limit E^g -^ E oi Eq. ^ and of the 
short-range approximation of Eq. (j4]) one can infer that 

^l;(0) = ^^^ f'(E) - With this the short-range approxi- 
mation of Eq. ([3]) can be written in terms of f{E) as 
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IE flho {f{E) - flbg/aho)^ 



7bg ahgfiAB 



f'iE) 
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Our model can be easily extended to anisotropic har- 
monic traps. The eigenenergy relation in a trap with 
ujx — ^y ~ rj'^z is known to be a = —^/TTd/J-{x,rj) 
with d, x, and J^ defined in [y]. Since a{E,B) is trap- 
independent we have to replace for r; 7^ 1 in the eigenen- 
ergy relation a\iof{E) by —y/Trd/ F(x,ri). One can show 
that this necessitates the same replacement in the ex- 
pression for $|;(0) and accordingly in Eq. (|9]). 

For the realistic case of ^Li-^^Rb we have performed 
full numerical multi-channel (MC) calculations in order 
to obtain eigenenergies and channel admixtures for dif- 
ferent trap frequencies w (77 = 1) and magnetic fields 
B. From the £' — )• limit in free space we obtain 
flo = — 17.77a.u. and from the Cg coefficient _Rcff = 
1899.9 a.u. The magnetic field positions of vanishing and 
resonant scattering length and the channel admixtures 
at resonance of the first two trap states in a shallow trap 
with Lj — 2tt X 20 kHz yield the parameters /i = iAAfiB, 
a* = 63.14 a.u., Bq = 1064.62 G, je = 2.73 x lO^^^a.u., 
and 7bg = 2.63 x lO^'* a.u. 

Figure [1] shows a comparison of the eigenenergies and 
RBS admixture obtained from the full MC calculation 
and from our model in a trap with w = 27r x 200 kHz 
which corresponds to a deep optical lattice. Both results 
are in very good agreement with a deviation < 0.003 hu 
and < 0.1%, respectively. This shows that the model ac- 
curately covers the E and B dependence of the scatter- 
ing process. Only for energies well below zero, the model 



fails to reproduce E and A correctly. Here, the van-der- 
Waals interaction becomes dominant the long range such 
that the approximation of a short-rang interaction breaks 
down. For the considered trap the roles of a* and a^g [E) 
become apparent through a significant broadening of IS.B 
by 0.18 G between the 1*^* and the 4*^ state. 
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FIG. 1. a) Energy spectrum of ^Li-*^Rb as a function of the 
magnetic field B in a trap with to — 2-71 x 200 kHz. Dots indi- 
cate MC calculations while lines indicate solutions of Eq. ^ . 
The eigenenergies of Eq. (|lip are shown as dashed lines, b) 
Scattering length a{E,B)/aho ~ f{E) corresponding to an 
eigenenergy E of the system. At _E = huj{2n -\- 1/2) with 
n — 'H resonances appear, c) Contribution A^ of the RBS for 
each energy level as a function of the magnetic field B. Dots 
indicate MC calculations while lines indicate results of Eq. (|9]) 
with energies from solutions of Eq. © . 

As shown in Fig. [1] b) the resonances of the scattering 
length a{E, B) ex f(E) are located at E^"^} = fiw(2n-t-i). 

In an anisotropic trap the resonance energies E^cs are 
determined by the roots of J-{x,ri). Consequently, the 
magnetic resonance position changes according to Eq. ([7]) 
from the free-space position Bfi = Bq — SB to 



B^:}^Bo-SB + E^:^/f,. 



(10) 



The difference of the resonance position for each energy 
level n opens the exciting possibility to tune the mag- 
netic field to a resonance of a specific trap state which in 
turn enhances inelastic collisions depopulating this level. 
By successively adjusting the magnetic field at different 
resonance positions one might be able to engineer an 
ensemble in an excited state or cool the system to its 



relative-motion ground state. A good candidate for this 
approach would be an MFR of ^'^^Cs at 19.8 G where the 
small value of fiAB — 2tt x A kHz 17[ admits to address 
single levels in reasonably deep traps. 

Applying Eq. (ITUl) one is able to explain the disagree- 
ment of an experimentally observed MFR position of 
s'^Rb in a negligibly weak trap {B.^s = 9.09(1) G [18]) 
and a trap of frequency w^ = Stt x 33 kHz, ujy ^ uJz ~ 
27r X 42kHz (B^cs = 9.121(9) G 0). The energy de- 
pendence of SB in unknown. However, its impact is 
likely to be negligible. For Ce = 4660 a.u. ^ it holds 
d5B/dE\^^Q < 0.1 /i p f or both Obg ^ a* and Obg ^ a* 
while l//i = 0.5 /ij^^ [i7|. Hence, the resonance shift is 
approximately given by E^d/^, = 0.034 G which is in 
good agreement with the experimental results. 

Another effect of the trap concerns the RBS admixture. 
It is present for each energy level [see Fig. [T] c)] which 
suggests that RBS_ molecules can be created not only in 
the bound state |2(| but also in exited states, e.g. exited 
bands of an optical lattice. Indeed, these confinement- 
induced molecules have been experimentally observed by 
Syassen et al. [9|. By inducing Rabi oscillations between 
atoms and RBS molecules at a very narrow ®^Rb res- 
onance {fiAB = 27r X 2 kHz, abg = 100.8 a.u.) they 
could produce a large number of molecules in an opti- 
cal lattice with two atoms per site in the center. After a 
sudden change of the magnetic field they measured the 
number of unbound atoms featuring pronounced maxima 
and broad minima. The suppressed dissociation at the 
minima can be attributed to a strong RBS admixture of 
excited trap states. Supported by our MC calculations 
for ^Li-*'^Rb we assume je ~ 7bg in order to estimate 
the RBS admixture using Eq. ([9]). Figure [2] shows the 
atom number observed in [9| and the RBS admixture for 
excited eigenstates at different magnetic fields. Clearly, 
a large RBS admixture coincides with suppressed dis- 
sociation. Here, the RBS molecules survive as part of 
an eigenstate of the new Hamiltonian while for magnetic 
fields where the RBS admixtures are small the projec- 
tion of the RBS onto the eigenstates leads to a larger 
fraction of unbound atoms. The effects of the energy- 
dependence of Obg can be studied by an effective range 
approximation abg = (oq ^ — jk^ReS + \k'^Vcs)^^ ■ With 
ao — 100.8 a.u., i?eff determined from C^ = 4660 a.u. [1^ 
and Vcs — (1700 a. u.)'^ the positions of small RBS admix- 
ture are shifted towards those of maximal dissociation. 

For no external trap the resonance of the scattering 
length coincides with the maximal RBS admixture to 
the scattering wave function [ll| such that the infiu- 
ence of both effects can hardly be distinguished. In 
traps this rule can be strongly violated. Searching for 
roots of the derivative of Eq. © with respect to E 
one finds a maximal RBS admixture where E solves 
abg(S)/aho = I{E) - 2f'{Ef/f"{E). For large Obg and 
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higher lying solutions i?max of this equation ii^Aiax can be 
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FIG. 2. Atom number vs. magnetic field (dots) as measured 
in [3|. RBS admixture of eigenstates according to Eq. Q 
with 7_B = 7bg, a* ^ Obg, and Obg = ao (solid line) and for 
an energy-dependent flbg given in the text (dashed line). 



(n) 

far away from the resonance energies Eics ■ Translated to 
the magnetic field the offset between the resonance po- 
sition Bres and the position of maximal RBS admixture 
Bmix can even approach AB which is accompanied by a 
vanishing scattering length at Bi^i^l Hence, this offset 
should be significant in a Fermionic system such as ®Li 
with a large background scattering length. 

Bourdel et al. [10] performed an experiment with 
2N — 7 X 10'* ^Li atoms in two different hyperfine 
states in a trap with oj^; = 27r x 0.78 kHz and u)y ss 
a;^ sa 27r X 2.2 kHz [lO[. They found a local maximum 
of atom loss close to Br but a global one at an about 
—80 G shifted magnetic field. Atoms at the Fermi edge 
have a relative-motion energy equal to the Fermi en- 
ergy Ep — 2TTh{6NLUxU}yUJz)^/^ . For aQ_= — 1405a.u., 



AB\e=o = -300 G, Cg 



1393.4 a.u. 



17| and 77 = 3 
our model predicts a maximal RBS admixture — 80.8G 
shifted from the resonance. This agrees well with the 
maximum loss position which can be an indication that 
the RBS admixture enhances transitions to deeper bound 
states and thereby influences atom-loss processes. Note, 
that another qualitative explanation for the off-resonant 
loss has been given by Bourdel et al. [10[ . 

Finally, we derive Landau-Zener transition propabili- 
ties for each avoided crossing in the spectrum. Expanding 
f{E) « f{Ei,g)+f{Ei,g){E-Ebg) inEq. ® around some 
background energy i?bg yields the eigenenergy equation 



{E-Eb- fiSB){E - E^g) = ^AB 



/(^bg) 
/'(^bg) 



(11) 



which describes the avoided crossing of a molecular 
eigenstate to a background state with coupling strength 
6"^ = /iAB/(^bg)//'(£'bg). For the n'^ avoided crossing 
and Obg < flho we have l/f'{Ebg) w 2fia;("^^f j/^F 
[13| such that the diabatic transition probability is 
given as P„ = p-27rG„ ^^^^^ q^^ ^ \5^/{h^lB)\ w 

2a;ABabg("^i(^)/(v^aho^) ■ Of course, only for |^| < 
fvjj the Landau-Zener theory can give exact results while 



otherwise two coupled states offer only a quantitative ap- 
proximation. This can be judged from Fig. [1] a) for the 
first avoided crossing where |(5| ~ lAhoj. 

In conclusion, we developed an analytic model of atoms 
in isotropic and anisotropic harmonic traps experiencing 
a Feshbach resonance. The energy-dependent scattering 
length was determined and compared to a previous pa- 
rameterization. Consequences of the model including a 
resonance shift, molecules in excited trap states, and a 
maximal molecular admixture away from the resonance 
were studied. Our model is in agreement with full nu- 
merical calculations and experimental results. We expect 
the approach to be applicable for an analytic treatment 
of other Feshbach-type resonances in a quasicontinuum. 
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